Sampling constrained probability distributions using Spherical Augmentation

Statistical models with constrained probability distributions are abundant in machine learning. Some examples include regression models with norm constraints (e.g., Lasso), probit, many copula models, and latent Dirichlet allocation (LDA). Bayesian inference involving probability distributions confined to constrained domains could be quite challenging for commonly used sampling algorithms. In this paper, we propose a novel augmentation technique that handles a wide range of constraints by mapping the constrained domain to a sphere in the augmented space. By moving freely on the surface of this sphere, sampling algorithms handle constraints implicitly and generate proposals that remain within boundaries when mapped back to the original space. Our proposed method, called {Spherical Augmentation}, provides a mathematically natural and computationally efficient framework for sampling from constrained probability distributions. We show the advantages of our method over state-of-the-art sampling algorithms, such as exact Hamiltonian Monte Carlo, using several examples including truncated Gaussian distributions, Bayesian Lasso, Bayesian bridge regression, reconstruction of quantized stationary Gaussian process, and LDA for topic modeling.Comment: 41 pages, 13 figure

Coupling of alpha(1)-Adrenoceptors to ERK1/2 in the Human Prostate

Introduction: alpha(1)-Adrenoceptors are considered critical for the regulation of prostatic smooth muscle tone. However, previous studies suggested further alpha(1)-adrenoceptor functions besides contraction. Here, we investigated whether alpha(1)-adrenoceptors in the human prostate may activate extracellular signal-regulated kinases (ERK1/2). Methods: Prostate tissues from patients undergoing radical prostatectomy were stimulated in vitro. Activation of ERK1/2 was assessed by Western blot analysis. Expression of ERK1/2 was studied by immunohistochemistry. The effect of ERK1/2 inhibition by U0126 on phenylephrine-induced contraction was studied in organ-bath experiments. Results: Stimulation of human prostate tissue with noradrenaline (30 mu M) or phenylephrine (10 mu M) resulted in ERK activation. This was reflected by increased levels of phosphorylated ERK1/2. Expression of ERK1/2 in the prostate was observed in smooth muscle cells. Incubation of prostate tissue with U0126 (30 mu M) resulted in ERK1/2 inhibition. Dose-dependent phenylephrine-induced contraction of prostate tissue was not modulated by U0126. Conclusions: alpha(1)-Adrenoceptors in the human prostate are coupled to ERK1/2. This may partially explain previous observations suggesting a role of alpha(1)-adrenoceptors in the regulation of prostate growth. Copyright (C) 2011 S. Karger AG, Base

PeakRegressor Identifies Composite Sequence Motifs Responsible for STAT1 Binding Sites and Their Potential rSNPs

How to identify true transcription factor binding sites on the basis of sequence motif information (e.g., motif pattern, location, combination, etc.) is an important question in bioinformatics. We present “PeakRegressor,” a system that identifies binding motifs by combining DNA-sequence data and ChIP-Seq data. PeakRegressor uses L1-norm log linear regression in order to predict peak values from binding motif candidates. Our approach successfully predicts the peak values of STAT1 and RNA Polymerase II with correlation coefficients as high as 0.65 and 0.66, respectively. Using PeakRegressor, we could identify composite motifs for STAT1, as well as potential regulatory SNPs (rSNPs) involved in the regulation of transcription levels of neighboring genes. In addition, we show that among five regression methods, L1-norm log linear regression achieves the best performance with respect to binding motif identification, biological interpretability and computational efficiency

Significance testing in ridge regression for genetic data.

Published versio

Behavioral genetics and taste

This review focuses on behavioral genetic studies of sweet, umami, bitter and salt taste responses in mammals. Studies involving mouse inbred strain comparisons and genetic analyses, and their impact on elucidation of taste receptors and transduction mechanisms are discussed. Finally, the effect of genetic variation in taste responsiveness on complex traits such as drug intake is considered. Recent advances in development of genomic resources make behavioral genetics a powerful approach for understanding mechanisms of taste

BALB/c Mice Deficient in CD4+ T Cell IL-4Rα Expression Control Leishmania mexicana Load although Female but Not Male Mice Develop a Healer Phenotype

Immunologically intact BALB/c mice infected with Leishmania mexicana develop non-healing progressively growing lesions associated with a biased Th2 response while similarly infected IL-4Rα-deficient mice fail to develop lesions and develop a robust Th1 response. In order to determine the functional target(s) for IL-4/IL-13 inducing non-healing disease, the course of L. mexicana infection was monitored in mice lacking IL-4Rα expression in specific cellular compartments. A deficiency of IL-4Rα expression on macrophages/neutrophils (in LysMcreIL-4Rα−/lox animals) had minimal effect on the outcome of L. mexicana infection compared with control (IL-4Rα−/flox) mice. In contrast, CD4+ T cell specific (LckcreIL-4Rα−/lox) IL-4Rα−/− mice infected with L. mexicana developed small lesions, which subsequently healed in female mice, but persisted in adult male mice. While a strong Th1 response was manifest in both male and female CD4+ T cell specific IL-4Rα−/− mice infected with L. mexicana, induction of IL-4 was manifest in males but not females, independently of CD4+ T cell IL-4 responsiveness. Similar results were obtained using pan-T cell specific (iLckcreIL-4Rα−/lox) IL-4Rα−/− mice. Collectively these data demonstrate that upon infection with L. mexicana, initial lesion growth in BALB/c mice is dependent on non-T cell population(s) responsive to IL-4/IL-13 while progressive infection is dependent on CD4+ T cells responsive to IL-4

Epidemic History and Evolutionary Dynamics of Hepatitis B Virus Infection in Two Remote Communities in Rural Nigeria

BACKGROUND: In Nigeria, hepatitis B virus (HBV) infection has reached hyperendemic levels and its nature and origin have been described as a puzzle. In this study, we investigated the molecular epidemiology and epidemic history of HBV infection in two semi-isolated rural communities in North/Central Nigeria. It was expected that only a few, if any, HBV strains could have been introduced and effectively transmitted among these residents, reflecting limited contacts of these communities with the general population in the country. METHODS AND FINDINGS: Despite remoteness and isolation, approximately 11% of the entire population in these communities was HBV-DNA seropositive. Analyses of the S-gene sequences obtained from 55 HBV-seropositive individuals showed the circulation of 37 distinct HBV variants. These HBV isolates belong predominantly to genotype E (HBV/E) (n=53, 96.4%), with only 2 classified as sub-genotype A3 (HBV/A3). Phylogenetic analysis showed extensive intermixing between HBV/E variants identified in these communities and different countries in Africa. Quasispecies analysis of 22 HBV/E strains using end-point limiting-dilution real-time PCR, sequencing and median joining networks showed extensive intra-host heterogeneity and inter-host variant sharing. To investigate events that resulted in such remarkable HBV/E diversity, HBV full-size genome sequences were obtained from 47 HBV/E infected persons and P gene was subjected to Bayesian coalescent analysis. The time to the most recent common ancestor (tMRCA) for these HBV/E variants was estimated to be year 1952 (95% highest posterior density (95% HPD): 1927-1970). Using additional HBV/E sequences from other African countries, the tMRCA was estimated to be year 1948 (95% HPD: 1924-1966), indicating that HBV/E in these remote communities has a similar time of origin with multiple HBV/E variants broadly circulating in West/Central Africa. Phylogenetic analysis and statistical neutrality tests suggested rapid HBV/E population expansion. Additionally, skyline plot analysis showed an increase in the size of the HBV/E-infected population over the last approximately 30-40 years. CONCLUSIONS: Our data suggest a massive introduction and relatively recent HBV/E expansion in the human population in Africa. Collectively, these data show a significant shift in the HBV/E epidemic dynamics in Africa over the last century

Lifespan development of stimulus-response conflict cost: similarities and differences between maturation and senescence

Age gradient of the mechanism of stimulus-response conflict cost was investigated in a population-based representative sample of 291 individuals, covering the age range from 6 to 89 years. Stimulus-response conflict cost, indicated by the amount of additional processing time required when there is a conflict between stimulus and response options, follows a U-shaped function across the lifespan. Lifespan age gradient of conflict cost parallels closely those of processing fluctuation and fluid intelligence. Individuals at both ends of the lifespan displayed a greater amount of processing fluctuation and at the same time a larger amount of conflict cost and a lower level of fluid intelligence. After controlling for chronological age and baseline processing speed, conflict cost continues to correlate significantly with fluid intelligence in adulthood and old age and with processing fluctuation in old age. The relation between processing fluctuation and conflict cost in old age lends further support for the neuromodulation of neuronal noise theory of cognitive aging as well as for theories of dopaminergic modulation of conflict monitoring

Functional Interaction between Type III-Secreted Protein IncA of Chlamydophila psittaci and Human G3BP1

Chlamydophila (Cp.) psittaci, the causative agent of psittacosis in birds and humans, is the most important zoonotic pathogen of the family Chlamydiaceae. These obligate intracellular bacteria are distinguished by a unique biphasic developmental cycle, which includes proliferation in a membrane-bound compartment termed inclusion. All Chlamydiaceae spp. possess a coding capacity for core components of a Type III secretion apparatus, which mediates specific delivery of anti-host effector proteins either into the chlamydial inclusion membrane or into the cytoplasm of target eukaryotic cells. Here we describe the interaction between Type III-secreted protein IncA of Cp. psittaci and host protein G3BP1 in a yeast two-hybrid system. In GST-pull down and co-immunoprecipitation experiments both in vitro and in vivo interaction between full-length IncA and G3BP1 were shown. Using fluorescence microscopy, the localization of G3BP1 near the inclusion membrane of Cp. psittaci-infected Hep-2 cells was demonstrated. Notably, infection of Hep-2 cells with Cp. psittaci and overexpression of IncA in HEK293 cells led to a decrease in c-Myc protein concentration. This effect could be ascribed to the interaction between IncA and G3BP1 since overexpression of an IncA mutant construct disabled to interact with G3BP1 failed to reduce c-Myc concentration. We hypothesize that lowering the host cell c-Myc protein concentration may be part of a strategy employed by Cp. psittaci to avoid apoptosis and scale down host cell proliferation